Energy harvesting by utilization of nanohelices 
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An energy harvesting device based on nanohelices is presented. The energy harvesting scheme 
based on nanohelices involves the same rectification circuitry found in many household electronic 
goods, which converts alternating current (AC) from a wall outlet into a direct current (DC) supply. 
The presented device, however, involves the rectification of ambient electromagnetic waves rather 
than the AC source from a household wall outlet. 



PACS numbers: 41.20.-q, 42.1 



INTRODUCTION 



Antennas are common in many wireless devices, such 
as cordless phones, radios, and television sets. For ra- 
dio and telecommunication applications, antennas are de- 
signed to receive electromagnetic waves in the gigahertz 
frequency ranges. The electromagnetic signals received 
are converted into electrical currents, which in turn gen- 
erate sound, images, and so on depending on the type of 
device. In principle, antennas are the most fundamental 
energy harvesting devices. 

The idea of collecting solar energy by antenna dates 
back as early as 1970's(author?) jlj. Since then, re- 
searches in solar energy collection by antennas have 
slowly progressed(author?) 0, 0, S|- However, due 
to the limitations on physical size of antenna, it was only 
recently a significant achievement in energy harvesting by 
antennas has been realized for the infrared (IR) spectrum 
of electromagnetic waves (author?) For telecommu- 
nication applications, dimension of antenna is on the or- 
der of centimeters. For the IR spectrum of electromag- 
netic energies, the antenna size scales on the order of 
sub-microns and this makes harvesting energy from light 
by antenna even more challenging. 

The efficiency of an antenna strongly depends on its 
size (author?) 0. With the advent of nanotechnology, 
the abundance of sub-micron sized structures which can 
be used as antenna exists today. Nano structures, such 
as nanorods, nanotubes, and nanodots, are beginning to 
shed some light on harvesting energy from electromag- 
netic radiation in IR to ultraviolet (UV) spectrum of 
range (author?) 0, i, [lo| . The size (or dimension) is 
not the only physical property that affects the efficiency 
of an antenna. For more sophisticated antennas, its geo- 
metrical configuration, e.g., shape, significantly affects 
the antenna efficiency (author?) [7||. Antennas based 
on simple nanorods, nanotubes, or nanodots leave little 
room for manipulating their geometrical configurations 
for optimizations. This put heHcal antennas out of the 
picture for developing solar cells based on antenna the- 



ory. This is all about to ch ang e with recent developments 
in nanohelices(author?) [urH E, Q H [3 ■ 

The helical antennas are widely deployed technol- 
ogy, which is well documented and studied in litera- 
ture(author?) (itI . [T5 |. Perhaps, the most widely de- 
ployed helical antenna, but which is also least likely to 
be thought of as one, i.e., as an helical antenna, is the 
transformer found in many electronic appliances. Be- 
cause majority of battery operated electronic devices run 
on direct current (DC) power, the electromagnetic en- 
ergy harvested by an antenna, which is an alternating 
current (AC) power, must be rectified. The circuitry 
that rectifies an AC into a DC power is referred to as 
a rectifier and one of the breed of rectifiers, a half- wave 
rectifier, is illustrated in Fig. [TJ 
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Figure 1: Half- wave rectifier with capacitor. 



The AC source in the rectifier circuit is physically rep- 
resented by a transformer. Fig. [2l The primary and 
secondary coils are labeled as Pc and Sc, respectively. 
Basically, transformer is an union of two helical conduc- 
tors in close proximity, where one of them is termed pri- 
mary coil and the other is termed secondary coil. Making 
analogy with the radio station and a radio, the primary 
coil plays the role of radio station and the secondary coil 
plays the role of a radio. In a transformer, primary coil 
transmits electromagnetic waves and the secondary coil 
receives those waves and converts them into electrical 
currents. 
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Figure 2: Representation of a transformer as the AC source. 

The term rectification becomes meaningless unless the 
output voltage, Moad, and the input voltage, T4n, satisfy 
the rectification condition, 



l^n > V^. 



load 



Vv 



be; 



(1) 



where Vte is the voltage drop across the diode and VJn 
is the voltage amplitude of the AC source[26|. The 
smaller magnitude for Vbe translates into higher sys- 
tem efficiency. For normal diodes, Vbe ranges between 
0.6 ~ 1.8 volts. Schottky diode, which is a special type 
of diode with very low forward- voltage drop, has the Vbe 
between approximately 0.1 ~ 0.5 volts. 

The household wall outlet supplies anywhere from 
115 V to 220 V for VJn, which is much larger than the 
diode forward voltage drop. Therefore, Vbe of 0.6 ~ 1.8 
volts is not of much concern when rectifying power line 
voltage. But, is this the same case for rectification in- 
volving ambient electromagnetic waves? To answer this, 
I shall consider a simple rectenna illustrated in Fig. [3l 
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For the nanorod of length Z, the electromotive force 
(EMF) generated inside of it would be given by 



EMF 



= / Ep • dl = Epl, 



where, for simplicity, Ep has been assumed to be parallel 
to the length of nanorod. For the nanorod of 1 um in 
length, intensity of 1 W m~^ generates 



EMF 



2.7 X 10"^ V 



inside the nanorod. The Vemf generated inside the 
nanorod rectenna can be identified with the VSn of Fig. [T] 
and the rectification condition, Eq. ([1]), gives 

2.7X 10-5V> l1oad + He- 

But, this cannot be satisfied for any Moad > even with a 
Schottky diode, which is known to have very low forward- 
voltage drop. Can Vemf be amplified so that the rectifi- 
cation condition, Eq. Jl]), is satisfied for sufficiently large 
Moad? The answer is yes and this involves the secondary 
radiation process, which is illustrated in Fig. [H 
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Figure 4: Amplification of Vbmf by secondary radiation pro- 



Figure 3: Simple nanorod rectenna, i.e., "rectifying antenna." 

An ambient electromagnetic plane wave with E-p = 
|jEp|| for the magnitude of its electric field part produces 
intensity / given by 

1 2 / cw3x 10*ms"\ 
I - -ceoiip, I ^ g 9 ^ g4 a2 jn-3 kg-\ 

where c is the speed of light in vacuum, Eq is the permit- 
tivity of free space, and Ep is the electric field(author?) 
[l9| . In general, the radiation from sun or nearby heat 
source does not form plane waves. However, if the longest 
dimension of the nanorod is comparable to the wave- 
length of incidence wave, the incidence wave can be ap- 
proximated as a plane wave(author?) The inten- 
sity of one watt per squared meters corresponds to the 



The irradiance from ambient source, which is indicated 
by Ep (electric field part) in the figure, induces radiation 
from the primary coil, Pc, which is indicated by Erad 
(electric field part only). Solution obtained by solving 
Maxwell equations shows that E^ad ^ Ep when the sec- 
ondary coil, Sc, is placed very close to the primary coil. 
Since Erad acts as the incidence wave for Sc, the Vemf 
is amplified in Sc by factor of Erad/Ep. 

The amplification of Vemf by secondary radiation pro- 
cess can be qualitatively understood by recalling the 
multiple-slit experiment with coherent fight source. Fig. 
m When plane waves pass through a multiple-slit plate, 
at distance x away from the slits, wavelets couple either 
constructively or destructively depending on the location 
of y and this results in bright and dark intensity patterns 
on screen. The cross-section of nanohelix, which forms 
the primary and the secondary coils in Fig. IH resembles 
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the multiple-slit plate (except here, the slit pattern is in 
ordered zigzag form) . 




Figure 5: Multiple-slit diffraction. 



That being said, the plane wave condition and the 
coherence of incidence wave is crucial to the amplifica- 
tion of Vemf by secondary radiation process. The out- 
door sun-light or the irradiance from light-bulb are not 
plane waves if plane waves are thought of as wave front 
with definite degree of coherence, which can be measured 
by the visibility of interference(author?) [2^. To put 
it simply, the degree of coherence is a measure of how 
perfectly the waves can cancel due to destructive inter- 
ference (or the opposite, measure of how perfectly the 
waves can add up due to constructive interference). The 
coherence was originally introduced in connection with 
Young's double-slit experiment in optics, where the in- 
terference becomes visible when Hght is allowed to pass 
through small aperture such as pin-hole and the effect 
becomes more pronounced with smaller pin-holes regard- 
less of the light source. Young's experiment justifies the 
use of plane wave input for nanohelices considered here 
as its height and winding pitch scales on the order of 
wavelength(27|. 



II. NANOTRANSFORMER ENERGY 
HARVESTING DEVICE 

A. Device structure 

The physical layout of energy harvesting device based 
on nanohelices (or nanotransformer (author?) [2ll . [2^ ) 
is illustrated in Fig. [6] and its equivalent circuit diagram 
is provided in Fig. [7l Borrowing the terminology from 
display technology, I shall refer to each element in diode 
layer (indicated by N-type and P-type square pairs in 
Fig. [H) and nanoheHces making contact with the diode 
element as a pixel. For clarity, the pixel is indicated by 
dotted rectangular three dimensional cube in Fig. [8l 
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Figure 6: Nanotransformer energy harvesting device layout. 



Equivalent circuit 




Incident EM •wave 



Figure 7: Equivalent circuitry. 



The proposed energy harvesting device based on 
nanohelices involves three major device layers (capacitor, 
diode, and nanohelix layers), which can be processed in- 
dependently and sandwiched together for the final prod- 
uct. Fig. [HI Such process allows manufacturing of the 
proposed energy harvesting device at large scales. 
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Three major layers of device 

^^^^^^^^^■B } Capacitor layer 

} Diode layer 

}Nanohelix layer 




Figure 8: Three major layers of nanotransformer energy har- 
vesting device. 

The diode layer can be first prepared as a one large 
sheet of a diode, which then can be patterned to form 
pixels of diodes. The nanohelix layer in current device 
scheme can be prepared by growing helices in normal 
direction of the substrate as illustrated in Fig. [9l Equiv- 
alently, nanohelices can also be spin coated directly onto 
the surface of a substrate. In this case, nanohelices would 
most likely be positioned with its length parallel to the 
surface of a substrate. To prevent diode layer from col- 
lapsing onto the substrate of nanohelix layer, transparent 
dielectric spacer such as Si02 nanoparticles may be dis- 
tributed on the substrate surface as illustrated in Figs. [9] 
andm 
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Figure 9: Schematic of grown nanohelices. 
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Figure 10: The cross-sectional view of nanohelix layer with 
diode layer placed over it. 



B. Operation principle 

For a successful rectification, the Vemf generated 
across each pixel must be amplified large enough to sat- 
isfy the rectification condition, Eq. (U). As already dis- 
cussed in previous sections, the amplification of Vemf is 
achieved thru the process of secondary amplification, re- 
call Fig. m The role of primary coil, i.e., Pc in Fig. 
m in current device scheme. Figs. [6] and [7l is played 
out by nanohelices belonging to the neighboring pix- 
els. Due to constructive and destructive interferences 
of wavelets originating from different nanohelices, there 
would be pixels receiving amplified radiation fields and 
there would be those pixels receiving virtually no radi- 
ation fields. Only those pixels positioned in locations 
where constructive interference occurs would generate 
large enough Vemf to meet the rectification condition 
and, eventually, participate in energy harvesting. 

The number of pixels, where a pixel is indicated by dot- 
ted rectangular three dimensional cube in Fig. [HI plays 
the key role in the energy harvesting based on nanohe- 
lices. As an illustration, assume that each pixel has a 
dimension of 10 um x 10 um for its surface area. In an 
ideal close packing, about 100 million such pixels would 
be able to fit in an area of 10 cm x 10 cm. Weber et 
al. have experimentally shown that ZnO nanowire can 
carry up to roughly 330 nA of electrical current before it 
snaps (author?) [23|. The current limiting resistors in 
Figs. SllGl and [7] prevents the overloading of nanohelices 
with too much current, thereby saving it from a break 
down. If assumed that each pixel contains single nanohe- 
lix and that each nanohelix carries electrical current of 
130 nA, this amounts to a total of 13 A out of the device. 
Of course, only those pixels positioned in locations where 
constructive interference occurs to generate Vemf large 
enough to meet the rectification condition contribute to 
the total current. But, even if one assumes that only 10% 
of 100 million pixels contribute in energy harvesting, this 
still yields the total of 1.3 A from the device. 

Having said enough about the potential of energy har- 
vesting based on nanohelices, the validity of working prin- 
ciples behind the proposed device depend heavily on the 
amplification of ambient electromagnetic fields by sec- 
ondary radiation process, thru which process generates 
Vemf large enough to satisfy the rectification condition. 
Since the amplification of Vemf depends on both mag- 
netic induction ( Bj-ad ) and electric field ( E^d ) parts 
of the electromagnetic radiation from the primary helix. 



Vemf « 



EradI 
Brad I 



one needs to quantitatively show that indeed Ej-ad and 
Brad get amplified significantly. 

Brad 1 1 Brad' 



rad 



IE, 
llEp 



ad I 



» 1, 



Bp 



Br 



> 1. 



(2) 



The quantitative verification of Eq. ^ involves the solv- 
ing of Maxwell equations and this is the task which I set 
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out to do in the next sections. 



III. THEORY 

All phenomena involving interaction with electromag- 
netic waves involve Maxwell equations. To keep the topic 
presented here self-contained, I shall briefly summarize 
the kind of manipulations and approximations assumed 
in obtaining the vector potential partial differential equa- 
tion (PDE), which marks the starting point for the rest 
of analysis throughout this presentation. 



A. Maxwell equations 

Maxwell equations, in the form independent of partic- 
ular system of units, may be expressed as 



V • E = ingp, V • B = 
dB 

V X E = -r;— , V x B 
at 



0, 

9 dt 



where p is the total charge density, J is the total current 
density, B is the magnetic induction, E is the electric 
fleld, and the positive constants g, 7, and r] depend on the 
particular system of units being adopted (author?) [23 |. 
If one assumes that both charge and current densities are 
specified throughout space and assume that E, B, p, and 
J vary in time as exp {iLut) , Maxwell equations may be 
re-expressed in an alternate form as 



V • E = ingp, V • B = 0, 
V X E = -iuriB, V x B = AttjJ - 
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E, 



where ui is the angular frequency. For a non-static case, 
where 7^ 0, the electric divergence relation, 

V • E = ingp, 

becomes redundant [2HI and the problem of electrodynam- 
ics is reduced to solving the following set of Maxwell 
equations in harmonic frequency domain. 



V • B = 0, 

V X E = —iLurjH, 



V X B 



47r7J + — E. 
9 



(3) 
(4) 

(5) 



automatically satisfies Eq. ^ and such vector A is given 
a name "vector potential." Substitution of Eq. ^ in Eq. 
(0]) gives 

V X (E + iLUTjA) = 0. 

The fundamental theorem of vector analysis tells us that 
any scalar field $ satisfies the condition V x V<i> = 0. 
And, this implies 

E + iuji]A = V$, 

where the sign of <& is arbitrary. However, because it has 
already been defined in literature that E = — V$ for the 
static limit, where w = 0, one chooses $ — $ for the 
scalar field and the previous relation becomes 



E = -iujrjA - V$. 



(7) 



Equation |(7]) automatically becomes the static limit ex- 
pression in the limit oj goes to zero. 

The concept of vector and scalar fields simplify electro- 
magnetic problem to solving of a single Maxwell equation 
(HI). Insertion of Eqs. ^ and ^ into Eq. JH) gives 



V X V X A = 47r7J 



lUJ"f 
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{-iuji]A - V$) 



Application of the vector identity, 

V X V X A = -V^A + V (V • A) 
transforms the previous relation as 



V^A - V (V • A) = -47r7J 



9 9 

After some rearrangements, I arrive at the expression. 



V V- A 



$ = V^A 
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A + 47r7J. (8) 



Since any A and $ satisfying Eq. ^ and Eq. ((7| , respec- 
tively, solves Eq. ^ , one is free to choose any convenient 
A and $ so that Eq. ([5]) becomes solvable. Choosing the 
following expression for $, 



$=^V.A, 

0^7 



(9) 



makes the left side of Eq. ^ to vanish. And, the electric 
field, utilizing Eq. can be expressed as 



B. Vector potential 

I proceed by seeking a vector field solution that simul- 
taneously satisfies Maxwell equations ([3]) thru |[5]). Any 
vector field A, which satisfies the condition 



B = V X A, 



(6) 



E 



V (V • A) 



(10) 



With <I> defined in Eq. the vector potential satisfies 
the following partial differential equation (PDE), 



V^A + g^A = -47r7J, 



K ^ujj— > 0, 



(11) 
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where the constant K has the physical impUcation of 
being the wave number. Equation ifTTI) is the well known 
Helmholtz equation and its solution is given by 



A(R) 



^ 7 



Jexp(-ii4:||R- Rg 
IIR-Rsll 



(12) 



where Rg is the position of current density, J = J (Rg) , 
and the volume integration is performed over entire re- 
gion containing the current source. 



IV. ANALYSIS 

A. Nanohelix 

The simplest solenoid is given by a non-planar helical 
curve depicted in Fig. [iTJ If ei, e2, and es denote a right- 
handed system of mutually perpendicular unit vectors, 
then every spatial points on filamentary finite helix can 
be represented by 



Rs 



Ro + R', 



(13) 



where Rq is the position vector defining the local origin 
O' and R' is the vector defining the position of current 
source relative to O' . In Cartesian coordinates, R' and 
Ro are given by 



R' 

Rn 



E 



Xr.. 



a sm s, 



= z' = 6s. 



(14) 



For the rest of the presentation, I shall designate the trio 
(ei, e2, 63) with (x, y, z) in respective order. Similarly, for 
the O' coordinates, I shall designate the trio (61,62,63) 
with {x',y',z') in respective order. The coordinates x' 
and y' describe a circle of radius a, and the z', which 
coordinate defines the height of finite helix, increases or 
decreases in direct proportion to the parameter s. The 
vertical distance between the coils, which is known as 
the pitch, equals the increase in z' as s jumps by 27r. The 
pitch is hence given by 



pitch = 27r6. 



(15) 



Assuming Rg is differentiable and that Rq does not 
depend on parameter s, the vector which is tangent to 
the curve defining the finite helix is given by 



dRg 
ds 



where 



-asms. 



T2 



4=1 



acoss. 



(16) 



(17) 



The total length, I, of filamentary finite helical curve is 
found by taking the line integration of associated tangent 
vector with respect to the parameter s. 



I 



T • ds 



ITII ds 



>-\/a 



-2a- 




Figure 11: Finite helix, where ei || x \\ x' , 62 || y \\ y' , and 
ea II z II 2'. The symbol || reads as "parallel to." 



where ds is the differential arc length of the finite he- 
lix segment. One thus obtains the upper limit for the 
parameter s, 



I 



< s < c. 



(18) 



With the parameter s, the entire height of finite helix is 
given by he or bl/^/a? -\- P as indicated in Fig. (fTTj) . 

The vector potential integral of Eq. I|12p is integrated 
over the entire volume containing the current sources. If 
the current sources are confined to a filamentary finite 
helix whose spatial curve is represented by Eq. (fT3|) . 
the helix may be partitioned into segments of finite but 
equal sizes as illustrated in Fig. [121 The filamentary 
wire forming finite helix has a length of I. Within the 
representation parametrized in s, each of N segments has 
length of c/7V, where c is defined in Eq. ifTS]) . The number 
of segments is arbitrary, as one can slice the wire into as 
many pieces as he or she wants to. To make the argument 
more concise, the filamentary wire is sliced into enough 
segments so that J for the segment is a constant and the 
entire segment is identified with its center, Rg„, on the 
finite helix. Then, for a detector placed at location R, 
the detected vector potential, which has been contributed 
from the segment n on the finite helix, gets approximated 
by the expression 



N 



S (Rg - Rg„) exp i-iK ||R - Rg 



-dV, 
(19) 

where the volume integration is over all space, the quan- 
tity c/N represents the length of segment, the constant 
a„ is related to the local current density J„ for the nth 
segment by J„ — cUn/N, and S (Rg — Rg„) is the Dirac 
delta function, which has the integral property. 



/(Y)5(Y-X)dF = /(X) 
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Application of the integral property of Dirac delta func- 
tion on Eq. (fTOl) yields the result 



A 



7Q!„cexp {--iK ||R — Rg. 
iV||R-R«„|| 



(20) 



The detector receives contributions from all N seg- 
ments of the finite helix, not just from the nth segment. 
Therefore, summing over the contributions from all N 
segments of the finite helix, I have 



N 



N 



n— 1 n— 1 



7a„cexp {~~iK ||R - Rs„H) 
7V||R-R.„|| 



Finally, in the limit the slices become finer and finer, it 
becomes 



JV 



= lim y 



7Q;„cexp {—iK ||R — Rg 



n=l 



iVllR- R,^ 



In the representation parametrized by s, one notices 
that c/N = As and, as N goes to infinity. As becomes 
infinitesimal, i.e.. As ds. Also, as the sHces get finer, 
what was the center point for the nth slice becomes the 
exact point for the slice, Rgn Rs- Similarly, what was 
an average current density within the slice becomes an 
exact current density for the point Rg, an ^ J (Rg) . 
Hence, 



N 



lim V/( 



N 



J f {x) ds, 



and the vector potential expression for the finite helix 
becomes 



7 



Jcxp(-iX||R- Rg 
IIR-Rsll 



-ds, 



(21) 



where Rg is defined in Eq. I|T3|) with its components given 
by Eq. As the current density, J, has not been de- 

fined, the vector potential integral in parametrized form, 
Eq. (|2T1) . cannot be evaluated. The current density 
source for the finite helix system is cast into a quanti- 
tative form in the next section. 



B. Induced current 

Consider an electromagnetic problem depicted in Fig. 
[T3j where a plane wave front is impinging on the finite he- 
lix whose configuration is describe by Eq. ifTS]). The real 
world solenoid, no matter how small, always has cross- 
sectional area of finite size which holds current respon- 
sible for induced electromagnetic radiation. Modeling a 
real world solenoid can be difficult due to the complica- 
tions arising from a finite thickness for the cross-sectional 
area. However, the mathematical modeling can be sub- 
stantially simplified by letting the cross-sectional area of 




Figure 12: Segmented finite helix. 
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Figure 13: Plane wave incident on finite helix. The Foynting 
vector is indicated by Sp. 



the wire to go to zero and, at the same time, letting the 
current density to go to infinity in such a manner that 
the fiux of current along the wire remains constant. For 
the current carrying helical wire modeled within such ap- 
proximation, Eq. ifTSj) suffices for the description of finite 
solenoid. 

The vector field A, whose solution satisfies the PDE 
of Eq. ifTTj) , arises as a result of induced current within 
the finite helix. This induced current inside a finite helix 
is due to the electric field component of impinging plane 
wave as illustrated in Fig. [I3l Assuming that wire form- 
ing finite helix can be represented by an isotropic ohmic 
conductor, the current density J is given by the Ohm's 
law. 



o- = fj (Rg) , 
Ep = Ep (Rg) 



(22) 



where a is the electrical conductivity, Ep is the polar- 
ization (electric field) of impinging plane wave front, and 
Rg is the position of current source. 

Electromagnetic waves have both electric and magnetic 
field parts, as illustrated in Fig. [141 For a nanohelix 



whose coil diameter is less than 50 nm, an IR range elec- 
tromagnetic plane wave passing through it would be per- 
ceived as a DC magnetic field switching between on and 
off modes at a rate of wave frequency. This comes about 
because the wavelength of electromagnetic waves in the 
IR range scales on the order of microns, which implies 
the vast majority of IR range electromagnetic waves are 
more than hundred times larger in their wavelength com- 
pared to the 50 nm diameter of the nanohelix. And, this 
argument strengthens with nanohelices of smaller diam- 
eter size (or with longer wavelengths, for example, mi- 
crowaves and radio waves). 
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Figure 14: The Ep and Bp are electric and magnetic parts of 
electromagnetic wave, and Sp is the Poynting vector. The en- 
closed box shows the time varying flux of magnetic induction 
enclosed by the nanohelix loop, where the symbol x denotes 
the case of Bp = — Bpea and the symbol • denotes the case 
of Bp = -Bpes. 



The macroscopic conducting coil driven by a house- 
hold line voltage alternating at 60 Hz generates elec- 
tromagnetic waves radiating at 60 Hz. Such radiation 
would have a wavelength of roughly 5000 km in air, i.e., 
A = 3 X 10^/60. A secondary conductive coil placed 
nearby it would perceive as if the magnetic part of the 
radiating electromagnetic wave was a DC magnetic field 
switching in an empty between on and off modes, except 
now at the rate of 60 Hz cycle. Again, the argument holds 
because even a solenoid of diameter as large as 100 cm is 
still fifty million times small compared to the wavelength 
of A = 5000 km. 

In summary, the case involving nanohelix and the IR 
spectrum of electromagnetic waves, in principle, is just 
the scaled down analogue of the case which involves the 
interaction between electromagnetic waves radiated by 
macroscopic conductive primary coil driven by a house- 
hold Hue voltage alternating at 60 Hz and the secondary 
solenoid placed nearby. 

That being said, under the exposure of electromag- 
netic radiation, net electric current gets induced inside 
the nanohelix in accordance with Faraday's law of in- 
duction and the total induced current density inside the 
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Figure 15: Perpendicular and parallel polarization compo- 
nents at the local helix segment. 



nanohelix must be expressed as 

J = crEp Jb, 

where Jb is the induced current density contribution aris- 
ing from the magnetic part of the incidence ambient elec- 
tromagnetic wave. Nevertheless, for the analysis here, I 
neglect Jb as this involves very lengthy derivation on its 
own. However, it is reminded that Jb only makes J big- 
ger. Therefore, once I show that Eq. ^ is satisfied even 
with contribution from Jb neglected in J, redoing the 
problem with contribution from Jb included in J should 
only make the case even firmer. 

Returning from a short detour, in the case where finite 
helix is made of a filamentary wire, only the component 
of electric field which is parallel to the local length of 
wire can induce current as depicted in Fig. [TH The 
polarization of impinging plane wave can be decomposed 
into Epj^ and Ep|| at the local point Rg of the finite helix, 
where Epj^ and Epy are the two components of Ep that 
are, respectively, perpendicular and parallel to the local 
tangent of finite helix at Rg. In a filamentary wire, only 
Ep|| can result in induced current. Mathematically, Epy 
at local point, Rg, is expressed as 



E. 



Ep-T 
T T 



E 



p — Ep (Rg) , 
T = T (Rg) , 



(23) 



where T is the local tangent vector for finite helix at Rg. 
The expHcit expression for T has already been defined 
in Eq. Insertion of Eq. ^ for Ep in Eq. ^ 

ensures that only Epy takes part in the generation of 
locally induced current density. 



Ep-T 
T T 



T, J = J(Rg) 



(24) 



Without loss of generality, the electric field of imping- 
ing plane wave at points on the finite helix can be ex- 
pressed as 



Ep (Rg) = exp (zK • Rg 



iujt) ^ E, 



(25) 



where the wave vector K is given by 



Eq. (|32|) may be simplified to 



K = ^ K,e, 



In terms of the direction cosines, 



cos ai 



El 



pi 



i = 1,2,3, 



(26) 



(27) 



ii-^pi 

the Eq. l(25|) can be expressed as 



Ep — Ep exp (iK • Rg + iut) cosaiGi. 

1=1 

With the following expressions, 

Ep • T = EpX exp {iK ■ Rg + iujt) , 

where 

X = — asinscosai + a cos s cos 012 + 6 cos as, 
the J of Eq. |[24]) becomes 

jTD ^ ^ y exp (iK • Rg + iujt) e, 



(28) 



(29) 



(30) 



i=l 



with the superscript TD denoting the time domain. In 
the frequency domain analysis Eq. ([30]) simplifies to be- 
come 

J™ = ^-jH^^r,exp(*K.R,)e„ (31) 



where the superscript FD now denotes the frequency do- 
main analysis, and Ti and x are respectively from Eqs. 
dm and l(29ll. 

It helps to simplify the analysis in the proceeding sec- 
tions by re-expressing Eq. (|29|) in an alternate form. It is 
well known in mathematics that any linear combination 
of sine waves of same period but different phase shifts is 
also a sine wave of same period, but different phase shift. 
It can be shown then 

oi sin ip + a2 sin {(p + bi) 



of + a2 + 2aia2 cos bi 



Since 



(f + arctan 



sin {(fi + bi) — cos (p, 



a2 sin 61 
fli + a2 cos bi 



5i =±(n+f) 
n = 0,l,2,--- 



(32) 



ai sin if + a2 cos f — y af + 



ip + arctan 



(33) 



for bi ^ tt/2. Equation l(29|) is compared with Eq. ([33] 
to yield 



X — CL\/ cos^ ai + cos^ a2 

' cos a2 



X sm 



s — arctan 



6 cos as, (34) 



^cos ai 

where the odd property of arc tangent function, 

arctan (—6*) = — arctan 0, 
has been utilized in the final step. 

C. Induced fields 

1. Induced vector potential 
In Cartesian coordinates, one has 



IIR-Rs 



^ ^ (Xj Xq^ x^) 



i=l 



1/2 



where R = ^^'^ is from Eq. lfT3|) . With 

Eq. ([Ml substituted in, u becomes 



(xi — a^oi — acoss)^ + {x2 — Xo2 — asins)^ 



+ (2^3 - 2^03 - bsY 



Insertion of Eqs. ^ and ^ into Eq. ^ yields 



(35) 



(36) 



where 



G^ = [ — exp {iK ■ Rs - iKu) ds. (37) 
Jo u 



The quantity K • Rg in Eq. (|37|) is given by 



K • Rs = ^ {K.xoi + Kix[) , 



where 



Kix'^ = aKi cos 3 + aK2 sin s + bK^ 



3S. 



(38) 



(39) 
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Equation l(39|) is compared with Eq. l(33|) to yield 

3 



Kix[ ^a^Kl+Kl sin 
and Eq. (|38l) hence may be expressed 



s + arctan 



+ a \J + K2 sin s + arctan (j^ 



With Eq. igo]), the of Eq. ^ becomes 



■ exp {iv) ds, 



(40) 



(41) 



where 



V = KiXoi + K2X02 + K3 {xo3 + bs) - Ku 



aJKf + sin 



s + arctan 



(42) 



(43) 



Eq. (|4T|) may be separated into the real and imaginary 
parts 



Utilizing Euler's formula, 

exp (iv) — cos V + i sin v, 



3?Gi 
3G, 



xTi ( cost; 



ds, 



(44) 



to yield 



A = gjA + i3A, 



3?A \ 



5RG, 

QA y " a2 + 52 I 3G. 

2 — 1 



(45) 



Using trigonometric identities. 



cos Q cos — 2 ["^^^ ~ '^) + cos (6* + ip)\ , (52) 

sin d sin = 2 I*^*-"^ ^ '*') ~ cos {6 + (^)] , (53) 

sin B cos V ~ ^ [sin (^^ ~ v) + sin (f? + (,c)] , (54) 

the xTi of Eq. gSl) and XT2 of Eq. gZl) may be re- 
expressed into a canonical form, 

xTi = 2'^^*^! fee's (2s — C2) — cos C2] — a6c3 sin s, (55) 

XT2 = 2 "^^1 [^^'^ (^'^ ~ C2) — sin C2] + a6c3 cos s, (56) 



where ci, C2, and C3 are defined in Eqs. (|49|) thru l(5T|) 
and the even and odd properties of cosine and sine, 

cos (—a;) — cos x, sin (— x) — — sin a;, 

have been utilized in the result. The xTs of Eq. (|46l) . 

of Eq. ((55|) , and of Eq. (|56l) are substituted into 
Eq. (im) to yield 



aci . ^ , 2aci cos C2 

cos (i; + 2s — C2 ) cos v 

u u 



5RGi = V 







H cos (ti — 2s + C2) sm iv + s) 

u u 

^2&C3 . . J , 

H sm (w — s) as. 



(57) 



5RG2 



4io 
26c3 



2&C3 2aci . 

cos {v + s) sm C2 cos t; 

u u 

aci 



cos (w — s) H sin {v + 2s — C2) 

u u 



aci . , ^ s 

sm [v — 2s + C2) 

u 



ds, 



(58) 



In explicit form, the xTi in Eq. (|44|) . with the aid of Eqs. 
(flTl) and l(34|) . for each z = 1, 2, and 3, becomes 

xTi = —c?c\ sin (s — C2) sin s — a6c3 sin s, (46) 
XT2 = a^ci sin (s — C2) cos s + a6c3 cos s, (47) 
xTs = a^ci sin (s - C2) + 6^C3, (48) 



5RG3- 







26c3 ac\ . , , 
cos V H sm (z; + s — C2) 

u u 



aci 

sm (w — s + C2) 

u 



ds, 



(59) 



where, 



ci = \/ cos^ ai + cos^ a2, 
'cos a2 



C2 = arctan 



C3 = cos as. 



cos Qfi 



(49) 
(50) 
(51) 



3Gi-^ 







2fcc3 26c3 

cos (w + s) cos [V — s) 

u u 



aci 2aci 

sm (w + 2s — C2) cos C2 sm v 

u u 



aci . , „ , 

H sm (V — 2s + C2) 

u 



ds. 



(60) 
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3G2 



4io 



aci 



cos {v + 2s ~ C2) 



cos (w — 2s + C2 ) H sm (w + s) 

u u 

2aci . . 26c3 . , ,1 
sm C2 sm V H sm (v — s) as, (61) 







aci 

cos (w + s — C2) 

u 



aci bc3 . 

H cos ( u — s + C2 ) H sm v 

u u 



ds, 



(62) 



where ci, C2, and C3 are defined in Eq. (|49ll5ip . 

2. Induced magnetic induction 

Substitution of A = 5RA + i3A into Eq. © gives 
B V X 3fJ A + iV X 3A. 
Insertion of 5RA and 5A from Eq. l(45|) yields 

B = 3fiB + i!3B, (63) 

where 



5RB 



(64) 



curl 

The ith component of curl of Eq. ([64]) is given by 

3 



V X ^ f e, 



In terms of vector components, Eq. l(64|) becomes 

3 



3?B 
9B 



2—1 ^ ^ 



(65) 



where indices j and k are chosen in accordance with the 
cyclic rule. 



Ifi 



1 

2 I then j 
3 



3 I and k 

1 



1 1 , (66) 



and eijk is the Levi-Civita coefficient. 



+1, if {ijk) is (123) , (231) , or (312) 
e^jk = \ -1, if fefc) is (321) , (213) , or (132) 
0, otherwise 



(67) 



Expanding out the Levi-Civita coefficient e^fc, following 
the rule stated in Eq. the 5RB and 9B of Eq. feS]) 

become 



SRB \ jaEp ^ / a.sRGfe - dkftG, 
3B j a2 + 52 2^ I apCfe - dk^G, 



where the indices j and fc are assigned in accordance with 
the cyclic rule defined in Eq. I|66p . Equations (|57l) thru 
(|62| may be summarized in the following form. 



. Q 



Cct COS (w + 7?cO 



ds, 



(69) 

where ip represents SRGi or SGi, the sums J2g ^^'-^ 
denote summation over terms involving sines and cosines, 
respectively; and, (CsgXa) and (i?se,i^ct) are the respec- 
tive constant terms which can be identified from the ob- 
servation of sines and cosines of which involve v in their 
argument. The operator V in Eq. (|64l) operates only 
on the coordinates of the detector. Since only u and v 
involves the detector coordinates, it can be shown that 

5'V^ = /(E%V'«.-E^^")d., (70) 



where Z = a;, y, z or ^ = 1, 2, 3 and 



V'sg = COS (v + l?sg) div ■ 



4>u = sin {v + z?ct) div + 



cos {v + ??ct) diu 



To compute for diu, I utilize the vector operator identity, 
V||R-Rs||" =n||R-RsH""^(R-Rs). (71) 
With Eq. dn]), 5/ ||R - RJ" becomes 

nei ■ (R - Rs) 



a,|lR-RJ| =ei.V||R-RJ|" = 



R Rs 



,2-n 



or 



5/||R-Rs 



Since u = ||R — RgH , I have 



n{xi - Xsi) 
|R-R,II'^"' 



(72) 



Xl — X^l Ul 

oiu^ = , ui=Xsi-xi. (73) 

u u 



For the div, using Eq. l(42|) for v, one obtains 
Oiv = —Koiu = , 



(74) 



12 



where Eq. (|73|) has been substituted in for diu. Insertion 
of Eqs. ([TS]) and into Eq. ([70]) yields the expressions 
for di^Gi and 9/5Gi, 

where 



ma 



2 Jo 

aciK 



aciK 



cos {v + s — C2) + 



26c, 



■ cosw 



2hcaK 



cos (w — s + C2j H ^ sm (u + s — C2) 



sm V — sm (w — s + C2 J 



ds, (80) 



u/ifcos (w + T^sg) sin (u + i^sg) 



u/ cos (u + "i^ci ) u/-f^sin (u + t^ct) 



Comparing Eq. (|75l) with Eq. |(69]) , I obtain the following 
transformation rule, 



sin [v + dse) Ago, 

cos (w + l?ct) IpcLO, 



(76) 



where -dsg and 'i?ct are extracted from the argument of 
cosines and sines by direct comparison. The computa- 
tions of didiGi and 9;3Gi are done by simple replace- 
ments of sines and cosines in Eqs. (|57|) thru l(62|) following 
the rule defined in Eq. ((76l) . This yields the expressions 



dim. 



(77) 



where 



4 



cos (u -I- 2s — C2) ^ — cos (v -I- s) 



2aciC0SC2 2bcaK 

7: cos V -I 5 — cos [v — s) 

aci , ^ , aciK . , ^ ^ 
^5- cos (1; - 2s + C2) 7^- sm (w -I- 2s - C2) 

26c3 / 1\ / 2aciii'cosc2 
— — ( 1 H — ) sm yv + s) H sm v 



u 



-I 5- sm [v — s) — sm (u — 2s -|- C2) 



(78) 



5R*2 = T 



— ^ cos (u -I- 2s — C2) H 5- cos [v + s) 



2aci sm C2 2&C3 , , 
cos V -I — cos [V — s) 

ac\K aci . , „ , 
— 5— cos (V - 2s + C2) -I ^ sm (w -I- 2s - C2) 

2hcaK 2aciKsvac2 
— sm (w -I- s) H sm v 

2bc-,K 



— — sm [V ~ s) sm [v — 2s + C2) 



ds. 



(79) 



aciK 2bc3 . , 
— ^ cos (u + 2s — C2) -I 5- cos (v + s) 



2aciK COS C2 



2&C3 



■ cos V 



cos {v — s) 



aciK , „ , aci . , „ ^ 
— — cos (V - 2s + C2) H ^ sm (w + 2s - C2) 

2bc3K . 2aci cos C2 . 
— sm [v + s) sm v 



2bcaK , . , aci . , n \ 
^ T, — sm (w — s) H sm (u — 2s -t- C2) 



ds. 



(81) 



aci 

— ^ cos [V + 2s — 02) 



-I 5 — cos (V + s) 



2aciK siiiC2 



■ cosw 



2bcaK , . aci 
H ^ — cos (u — s) H ^ cos (w — 2s -I- C2) 

If' 

aciiiT . . „ ^ 2&C3 . , , 
-I 5— sm (w 4- 2s — C2) H ^ sm [v + s) 



2acisinc2 . 2fec3 . 



aciK 

V? 



smw -|- 



sin (i; — 2s -|- C2) 



■ sin {v — s) 



ds, 



(82) 



3 = 7: 







aci bc^K 
cos (v + S ~ C2) TT- 



aci 
w 
be. 



COS (W — S -I- C2 

aciK , 



aciK , 

— sm (V + s ~ C2) 



■ sm V 



sin (u — s + C2) 



ds, 



(83) 



where ci, C2, and C3 are defined in Eqs. (|49l) thru (|5T1) . 

The components of 5RB and SB are readily extracted 
from Eq. (|68l) to yield 



5RB 



E 

4=1 



where 



55, 



(84) 
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Utilizing Eq. the 'SiB, and QB, of Eq. ^ become 



a2 + 62 



(85) 

where « = 1,2,3 and the indices j and k are assigned 
in accordance with the cyclic rule defined in Eq. (|66|) . 
Knowing that Bp = Ep/co, where Co is the speed of light 
in vacuum, Eq. l(85|) may be expressed as 



^B, 
35, 



life 



(86) 

Since B = 3?B + z3B from Eq. ([Ml), the magnitude of 
B is given by 



B = VBtB = V(5RB - i^B) ■ (gfJB + iQ?B) 



= V3f?B • 5RB + 3?B • 3?B 



^ [(5Ri?,)' + {QB,y 



1/2 



With Eq. I|86p substituted in for ^Bi and one ob- 
tains 



7CrCo 



-Brad ^ _B_ _ 

Bp ^ + 62 



-Bp 



(wgSR^'l - Ul^R^'g)^ + (ui5R^'2 - U25R*l)^ 



(mi3*2 - U23*l) 



(87) 



where SR^'i and are from Eqs. thru ([831) and Ui 
is defined in Eq. (|73|) for each z = 1,2,3. 



where didiGi = didiGi since 9^/ {dxdy) = 9^/ (dydx) 
for any mixed partial derivatives (recall that notation 9^ 
or dx represents d/dx). Insertion of Eq. (|45)) into Eq. 
([891) finally yields 

E = KE + i3E, 

where 

3 3 

JRE = ujr^'^A + V V eAdi^Gu (90) 



QE = -wr/SRA - 



3 3 



^^e,9,az3?G,. (91) 



(a2 + 62) 

Utilizing Eq. for 9/30; and 9/3fiG/, one has 



didi ( jjj^' ) — di 



mi 



m. 



or 



diXi. 



Since the notation 9^ denotes d/dx for, say, « = 1, and 
notation xi denotes x for I = I, y for / = 2 and so on, 
one finds 

= otherwiL ■ (92) 

where 5^; is the Kronecker delta. The partial derivatives 
didi^Gi and didi^Gi hence become 



3. Induced electric field 

The associated electric field may be obtained from Eq. 
Wl- Insertion of Eq. ([36]) into Eq. ^ gives 



E = —iujTjA 



U) {c? + 62) 



The term involving Cartesian gradient operator V, where 
V = Y^i=\ ^i^i, can be expressed as 

(3 \ 3 3 3 

W-Y^eiGi] = V^a^Gi =^e,a,^9iGi 
1=1 J 1=1 i=l 1 = 1 

3 3 

= ^^e,d,diGi 



and Eq. (l88l) becomes 



E = —iuji]A 



igaEp 



3 3 



.(a2 + 62)^^ 



^^e.a.azGi, (89) 



and the SRE of Eq. ^ and 3E of Eq. ^ get re- 
expressed as 



5RE 



(tE, 



62 ^ 



joj'qQGi 



1=1 



a2 -I- 62 ^ 

i=l 



jLur]^Gi 



I,! 



1=1 



(93) 



(94) 



where Eq. (|45|) has been substituted in for 3f?A and SJA. 
To compute for di^"^i and di^'i'i, one notes that SR^*/ 
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and 9^'; can be summarized in form as 



where 



T = 



E 



- e 

t;2 



sin(« + ^bj^^ J^„g(^^^cj 



cos (w + 1?' 



(95) 



where T represents or Q'^';, the sums 
and denote summation over terms involving sines 
and cosines divided by or u^; and, {CsgXsaXcuCc-;) 
and -(Jg^, t?J;^, t?J?^) are the respective constant terms 
which can be identified from the observation of sines and 
cosines of which involve v in their argument. The op- 
erator di only acts on non-source coordinates, of course. 
Since u and v are the only terms with non-source coor- 
dinates, one has 



\ /-a 'ya I \ ^ /-b 'Y'b 



where i = x, y, z or « = 1, 2, 3 and 

cos -I- -d^^ diV sin (u -|- d^g) diU^ 



(96) 



Uj K CO^ (v + 



sin (t; + i^^^) 



u 


2 




UiK COS (w + -d]^^] 


2ui sin (w - 




u 






3ui cos {v + 


sin [v - 
u 




2ui cos {v + d^^) 


UiK sin - 






u 





Comparing Eq. (|98|) with Eq. I|95p . one identifies the 
transformation rule for the sines and cosines given by 



sin {v + ^-J^T-,, 
cos(t; + ^^J^T^„, 



sm (y 
cos (u -h 



"serf ^ <To; 



(99) 



where 'd^g, and can be identified by observing 

appropriate cosines or sines in expressions for 3^4'; and 
Q'i'i of Eqs. dlH]) thru Application of Eq. ^ on 

Eqs. ^ thru {SS]) yields 



COS 



{v + -d^^) d^v sin {v + d^^) diU^ 



3A/ 



(100) 



sin {v + i^jjj 9iU COS {v + i?J;J 



where 



sin (w + diV COS (u -I- diU^ 



Utilizing Eq. l(72|) . it can be shown 

— 2 (x^ ^si) : — 3 {xi ^si) ^ 

or, since Xi — x^i = — Uj, 

diU^ — —2ui, diU^ = —3uiU, 



(97) 



where the dummy index I has been replaced by another 
dummy index i, of course. The expression for diV has 
already been defined in Eq. (|74|) . i.e., let i = I. With 
Eqs. {Ill) and {nzl), the expression for diT of Eq. l(96l) 
becomes 



,,2 



b 



E 



AC Ad ^ 



(98) 



3?Ai = 



aci 



2hc?.K ( + ) cos (u + s) 



cos (w + 2s — C2) 



-f- 2aci C0SC2 



3 



cosv 



+ 



+ 



if 2 3 

aci I r ) cos (v ~ 2s + C2) 



iaciK . , n N 
. ■ cos [v — s) -. — sm (v + 2s — C2) 



2 3 

2&C3 ( — - -4 - ) sin (w s) 



6aciifcosc2 . 



3 

sinu — 2bc3 ( — ^ ) sin (w — s) 



7 — sm (w - 2s + C2) 



ds. 



(101) 



15 



mo 



— cos (v + 2s- C2) 



2bc3 — cos (w + s) 

2aci sm C2 — ^ ? cos v 



^2 3 

2^C3 ( — - -5 ) COS - s) 



3aciK 

2 COS [v - 2s + C2) 



6bc:iK 



■ sm (V + s) 



•^2 3 
aci ( =■ ) sin (w + 2s — C2) 



GaciK smc2 . 



&bc^K 



sin {v — s) 



~aci 



sin [v — 2s + C2) 



ds. 



(102) 



a 

GbcsK 



aci 1 — =■ I cos {v + 2s — C2) 

6aciXsinc2 

. cos (u + s) 1 cos V 

4 ^,4 



aci — cos (u — 2s + C2) 

%bc-iK 3aciK . . , „ , 
■ -. — cos (u — sj H 1 — sm [V + 2s — C2) 

(K^ 3 \ . , 

2OC3 ^ p sm (v + s) 

\ U'^ J 

(K^ 3 \ . 

2aci sm C2 — 5 r sm v 

\u-^ vP ) 

[K^ 3 \ . , 



-. — sm fu — 2s + C2 j 



(105) 







3aciK 

2 COS (w + S - C2) 



3\ Saciif 
20C3 — cos u — COS (w - s + C2) 



aci 



+aci 



3 \ . ^ ^ ^ 6^g3^ . 
— sm ( u + s — C2 ) — sm v 



sin (v — s + C2) 



ds, 



(103) 



2 Jo 



K'^ 3 
aci I — - — ) cos (^; + s - C2) 



3\ , 

cosw — aci — cos (u — s + C2) 

If* \ u-' M"" y 

3aciJC , u 3 \ . 
^ — sm 1 1; + s — C2 ) — 0C3 ^5 ] smv 

,4 \ lf5 I 



3aciK . 

2 — sm (V - s + C2) 



ds, 



(106) 



where ci, C2, and C3 are defined in Eqs. (|49|) thru (|5T|) . 
Insertion of Eq. I|100p into Eqs. ([Ml) and ([94]) yields the 
expression given by 



5RE 
3E 



^2 , /,2 



+ 62 ^ V 

2—1 



where 



3Ai = 



iaciK 

-. — cos (u + 2s - C2) 



4 Jo 

3 \ / , ^ 6aciii'cosc2 

— 26c3 cos (u + s) cos?; 

3 \ , 

+ 2bc3 cos (u - s) 

\u-^ ) 

Zac\K dbc^K . 
H -. — cos (v - 2s + C2) — sm [v + s) 

(K^ 3 \ . , ^ 

— aci ^5 ^ sm (w + 2s — C2) 



i^2 

+ 2aci cos C2 ( — 5 F I sin v + 



66c3ivr . 



-^2 3 
-aci ( —3 g ) sin (u — 2s + C2) 



sin (u — s) 



ds, (104) 



-E 



3A, 



(107) 



Direct expansion of "REi and for each i = 1, 2, 3 to 
yields 



^Ei = jurj'iGi 
, 9 



- (u?3Al + -U2Ul9A2 + M3Ul9A3 - Q-^i) 



(108) 



lkE2 = Tc^ry^Ga 



— (uiM23'Ai + U29A2 + U3W23A3 — 3^*2 



(109) 
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(jj 



(110) 



and 



- (m?3?Ai + U2U15RA2 + U3U15RA3 - mil) , 



(111) 



'^£2 = jujr/^G2 



(uiU23?Ai + uj^A2 + U3U25RA3 - ^2 



(112) 



3i?3 = jujrj^Gs 



- (uiMsSRAi + U2U35RA2 + uMAs - ^3 



(113) 



The magnitude of E is given by 



E = VEtE = Vl^f^E - i3E) • (3fiE + i3E) 
= V3fiE • 3?E + 3E • Q?E 

or with Eq. (|107p substituted in for 3ffE and 3E, I obtain 



EraA E a 



a^ + P ^^^^ 



^ [i^E^f + {^E,y 



1/2 



(114) 



where IREi and '^Ei are defined in Eqs. ifTOSj) thru (fTTSl) . 

V. RESULT 

The fields are measured along the surface of cyHndrical 
shell illustrated in Fig. [I6l Relative to the O' frame of 
reference, an arbitrary point on the surface of cylindrical 
shell is given by 

R:, = R' + R^d> (115) 
where R' is from Eq. (fT3|) and is given by 

R' — a cos (s) ei + a sin (s) 62 + bse^ (US) 

with ei = x', 62 = y', and 63 = z'. In terms of cylindrical 
coordinates (p', 4>' ,z'^) , R^ can be expressed as 

Rj = p cos {(j)') ei + p sin [<p') 62 + 2^63 

and Eq. (|115p may be solved for R^^ to yield 

Rgj — (p' cos (j)' — a cos s) ei + (p' sin (j)' ^ a sin s) 62 



+ (4 - ^s)e3, 



(117) 



Z' = N-AZj— t 
Z'=(N-1)-AZ' 




Figure 16: Position of detector relative to the primed origin. 

where p' is a constant for a fixed cylindrical shell, 4>' 
sweeps from to 27r, and z' = z'^ ranges from —<x> to 00. 

Relative to the frame of reference O, Fig. [17] in which 
frame the unit bases x,y,z satisfy the condition x \\ x' , 
y II y', z II z', the locations R^ and R' are given by 

Rd = Ro + Rd' I^s = Ro + R', 

where Rg = Rq + R' is from Eq. (fT3l) . These relations 
can be combined with Eq. IjllSp to yield 

3 

R-d = R-sd + R' + Ro = R-sd + R' + ^ XoiBi. 



With R' inserted from Eq. I|116p and R^^ substituted in 
from Eq. (|117p . Rd becomes 

Rd (a;oi + p' cos <j)')ei + {xo2 + p' sin 0') 63 

+ (a;o3 + 4)e3. (118) 

The Rd in cylindrical coordinates can be expressed as 

Rd = pcos^ei + psin0e2 + Zde3, (119) 

where p is now the radius with respect to the unprimed 
reference frame O. Combining Eqs. (|118p and (|119p . I 
obtain 

{p cos (j) — Xqi — p' cos (/)') ei 
+ {p sin 4>- Xo2- p' sin (j)') G2 

+ (^d - 2:03 - 4) es = 0. (120) 



Equation l|120p can only be satisfied if and only if coeffi- 
cients of ei, 62, and 63 are independently zero, hence 

p cos (j) — Xoi — p' cos (f)' — 0, 
p sin (j) — Xo2 — p' sin (j>' — 0, 

Zd - Xo3 - 4 = 0- 



The third relation readily gives 

Zd = Xo3 + z'. 



(121) 
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Figure 17: Position of detector relative to the unprimed ori- 
gin. 



and the constant K in Eq. (fTTI) gets identified as 



where the free space electric permittivity Sq and the mag- 
netic permeabihty /io have the value given by 

eo ~ 8.854 X lO"^^ s^ A? m"^ kg-\ 
= 47r X 10"'^ mkgs"^ A"^. 

The total electric field and the total magnetic induc- 
tion part of the electromagnetic wave reaching the surface 
of an imaginary cylindrical shell, illustrated in Fig. [iTl 
can be summarized as 



E 



total 



= E' 



E 



radi 



B 



total 



= Bp + B 



radi 



and the first two relations rearranged to give 



psm( 
pcosi 



Xo2 



p' sin <p' , 
p' cos (f)' . 



From the ratio of the two, i.e., tanc 
obtain 



cost 



arctan 



+ P sin ( 



(122) 



^a;oi + p' cos(t)' , 

The p is found by combining the two relations, i.e., 
p sin + p cos 0, to get 

_ + + p' (sin 0' + cos 0') 
^ sin + cos 

The sin a + cos a can be combined utilizing Eq. l(33|) to 
yield 

sin a + cos a = \/2 sin \a + arctan (1)] 
= \/2sin [a + ^) 

and p becomes 



Xq\ "t~ Xq2 

+ p'v^sin I 



V2sin((/.+ f) 
Insertion of Eq. (|122p for (/) yields the result (2£ 

a:^oi+a;o2 + p'\/2sin(0' + f) 



\/2sir 



arctan 



I Zo2+P' ^ 
\ a;oi+p' c 



(123) 



where < 0' < 27r. With Eqs. (fm]) and (fT23l) . the sur- 
face of cylindrical shell illustrated in Fig. [T6]is completely 
defined relative to the reference frame of O. 

In MKS system of units, where length is measured in 
meters, mass is measured in kilograms, and time is mea- 
sured in seconds, the constants 5, ry, and 7 of Eqs. ([3]) 
thru ^ are identified as 



9 



V 



1 



Mo 

47r 



where Ep and Bp respectively represent the electric 
field and the magnetic induction parts of electromagnetic 
waves other than the radiation from nanohelix reaching 
the detector. The measure of field amplification, there- 
fore, may be expressed as 



total 



EI3 + E, 



ad I 



otal 



IB' +B, 



ad I 



E. 



IB. 



In many situations, the magnitudes of Ep and Bp are 
only minute changes from that of incidence wave, i.e.. 



Ep 



^Pll ~ llEpll and jiBpi, ~ nx^p, 
pressions can be approximated as 



IB' 



and the previous ex- 



total 



IE. 



rad I 



Ev 

-Btotal 



l|Ep|| 

Brad! 



Bv 



IIB. 



> 



> 



Ep + Eradj 

l|Ep|| 

||Bp + Brad! 

IIBpII 



(124) 
(125) 



A successful rectification of ambient electromagnetic 
wave requires the field amplification criterion defined in 
Eq. d). 



E, 



total 



»1, 



total 



»1, 



where the subscripts in the magnitudes of electric field 
and magnetic induction have been modified from rad to 
total. Observing Eqs. (|124p and (|125p . the requirement 
for a successful rectification of ambient electromagnetic 
wave is given by 



Ep 



»1, 



B 



rad 



Bp 



»1, 



which is the field amplification criterion defined in Eq. 

I am now ready to plot the results. For convenience, 
the origins of two reference frames, O and O', were made 
to coincide each other. This makes Xoi = Xo2 = 0, Zd = 
z'^, and Azd = Az^. Furthermore, it had been assumed 
that the helix winding started at = = and the 
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vacuum was assumed for the medium holding both the 
finite heHx and the propagating incidence and radiated 
electromagnetic waves. That being said, Eqs. l(87|) and 
(|114p are computed at the surface of cylindrical screen 
of radius /?', Eq. I|123p . using Simpson method coded 
in FORTRAN 90 for numerical integration (author?) 
[iB i and assuming the following input values, 

wavelength A ~ 555.016 nm, 



27r 

wave number vector: K (K, 0, 0) = Kyi = -r-x, 

A 



7r 7r \ 

polarization: Ep ( ai = — , a2 = 0, as = — j = £'py, 



Magnetic induction amplification factor at Selectetl z^- (/i=lieig]it) 




Figure 19: The ratio Brad/Sp, which is measured on the sur- 
face of cylindrical screen of radius p' — 273 nm and sliced at 
indicated z^. 



radius of helix: a = 40 nm, 
pitch of helix: pitch = 50 nm, 
radius of screen: p = 273 nm, 
increment along z axis: Az^ — 0.0976 nm, 
helix conductivity: tj = 5 x 10^ S, 
fully streched length of helix: I = 5 um, 

helix winding number: ~ 19.5. 

The results are illustrated in Figs. [18] thru [211 where the 
parameter h in Figs. [HI and [20] represents the height of 
nanohelix, i.e., h = be, as illustrated in Fig. [iTJ The 
results illustrated in Figs. [T9|and[2Tl therefore, respec- 
tively represent slices of Figs. [HI and [20] at helix height 
indicated by z^- 



3DPlotof LOG,„(B„yBp) 




(^'^(degree) 



Figure 18: The ratio -Brad/Bp measured on the surface of 
cylindrical screen of radius p' = 273 nm. 



3DPlotofLOG,„(£,^y£p) 




(|)'^(degree) 



Figure 20: The ratio _£Jrad/^p measured on the surface of 
cylindrical screen of radius = 273 nm. 



Ekctric field amplificaiion factor al Selected Zj. (/i-heighl) 
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Figure 21: The ratio E^^^/Ep, which is measured on the sur- 
face of cylindrical screen of radius p' = 273 nm and sliced at 
indicated z^. 

For the particular case where the input electromagnetic 
wave is specified by the wave number vector K = if x and 
Ep = ii'py, the resulting induced radiation from a single 
nanohelix is characterized by an electric field part which 
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is distorted in profile. The distortion in electric field pro- 
file can be attributed to the local geometrical configura- 
tion of the nanohelix, e.g., the winding pitch, etc. At dis- 
tances far away from nanohelix, such effects arising from 
the geometrical configuration of nanohelix should be av- 
eraged out, resulting in a radiated electric field profile 
which is symmetrical in shape. To confirm such behavior 
for the scattered fields at distances far away from scat- 
tering nanohelix, Eqs. I|87p and l|114p were re-computed 
with only the radius of cyHndrical screen changed to 
p' — 233.04 um, which is about thousand times larger 
than the previous value for the radius, p' = 273 nm. The 
results are shown in Figs. [221 and [23] for the ratios in- 
volving electric and magnetic fields, respectively. As ex- 
pected, the profiles for both magnetic induction and elec- 
tric field portray a symmetry about common axis. 



3DPlotofLOG,„(i?„yBp) 




^''(degree) 



Figure 22: The ratio -Brad/Bp measured on the surface of 
cylindrical screen of radius p' — 233.04 um. 



3DPlotof LOG|||(£,,y£p) 
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The proposed energy harvesting device based on 
nanohelices is actually an array of many pixels, where 
the term pixel has been adopted from display technol- 
ogy to denote a smaller energy harvesting device unit 
composing of single nanohelix (or few nanohelices). In 
such many-body (many-nanohelix) system, contributions 
arising from interaction with different nanohelices must 
also be taken into account for a total effect, which is a 
very challenging task. However, without having to go 
through such difficult calculations involving many-body 
effect, one can be sure of the existence of locations where 
constructive and destructive interference of wavelets oc- 
cur within the layer containing nanohelices. The "num- 
ber" wins in energy harvesting device based on nanohe- 
lices. Since there are biUions of nanoheHces in the system, 
even if only one hundredth of them actually participate in 
the energy harvesting, the total output power generated 
would be substantial. And, the demonstrated calcula- 
tion based on single nanohelix, supports the possibility of 
harvesting energy utilizing nanohelices. Any calculations 
involving many-nanohelices effect should be deferred for 
the optimization stage of the development process. 



VI. CONCLUDING REMARKS 

The presented energy harvesting device based on 
nanohelices, in all respects, can be thought of as a 
miniaturized version of rectifier circuits with transform- 
ers found in many electronic systems. The only differ- 
ence is that rectifier based on nanohelices rectify ambi- 
ent electromagnetic waves, whereas the conventional rec- 
tifiers rectify AC source from the household wall outlet. 
As with all rectifiers, the rectification condition defined 
in Eq. JT]) must be satisfied before the proposed device 
can actually convert ambient electromagnetic waves into 
a useful DC electrical power. The rectification condition 
can be satisfied if the condition defined in Eq. ^ can be 
met. 



In this work, by utilizing the secondary radiation process, 
I have explicitly shown that the condition imposed by Eq. 
([2]) becomes feasible with nanohelices. 
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